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We compare the azimuthal correlations arising from three and two hadron production in high 
energy proton-proton and nucleus-nucleus collisions at ^/sjvjv = 200 GeV, using the leading or- 
der matrix elements for two-to-three and two-to-two parton-processes in perturbative QCD. We 
first compute the two and three hadron production cross sections in mid-rapidity proton-proton 
collisions. Then we consider Au + Au collisions including parton energy loss using the modified 
fragmentation function approach. By examining the geometrical paths the hard partons follow 
through the medium, we show that the two away-side partons produced in two-to-three processes 
have in average a smaller and a greater path length than the average path length of the away-side 
parton in two-to-two processes. Therefore there is a large probability that in the former processes 
one of the particles escapes while the other gets absorbed. This effect leads to an enhancement 
in the azimuthal correlations of the two-to-three with respect to the two-to-two parton-processes 
when comparing to the same processes in proton-proton collisions since in average the particle with 
the shortest path length looses less energy with respect to the away side particle in two-to-two 
processes. We argue that this phenomenon may be responsible for the shape of the away-side in 
azimuthal correlations observed in mid-rapidity Au + Au collisions at RHIC. 



PACS numbers: 25.75.-q, 25.75.Gz 

Many interesting phenomena have been observed in 
high energy heavy ion collisions at the Relativistic Heavy- 
Ion Collider (RHIC) since it began operation almost a 
decade ago. The suppression of the single hadron trans- 
verse momentum spectra in Au + Au collisions, as com- 
pared to normalized p + p collisions is one of such phe- 
nomena. This is believed to result from the energy loss 
of fast partons traversing the medium and multiply scat- 
tering (both elasticly and inelasticly) from medium con- 
stituents P, Q ■ This picture of energy loss was further 
confirmed by azimuthal correlation studies where for cer- 
tain combinations of leading and associated particles a 
disappearance of the away-side peak is observed. This is 
interpreted as the absorption of the away side parton by 
the medium. Combined with an enhancement of the sin- 
gle hadron production in d + Au collisions as compared 
to normalized p + p collisions, these observations estab- 
lished the presence of a final state medium produced in 
Au + Au collisions at RHIC. Furthermore, the strong 
collectivity of the medium, as measured by the flow pa- 
rameter V2, among other observations, have led to the 
consensus that the medium is partonic in origin as well 
as strongly interacting 0, Q ■ 

In spite of the success of the above interpretation of 
RHIC data ||, there still remain important aspects of 
energy loss dynamics which are poorly understood. The 
picture has been further complicated by the recent ob- 
servation in azimuthal correlation studies in Au + Au 
collisions, where for the case when the magnitude of 



the momentum difference between leading and associated 
particles increases, either a double hump structure or a 
broadening of the away-side peak appears, while these 
structures are absent in p + p collisions at the same en- 
ergy 0, 0| ■ This observation gave rise to a large number 
of theoretical explanations which are based on the as- 
sumption that collective phenomena are at work in A 
+ A collisions in contrast to the case of p + p colli- 
sions @, |, G3, El Ei El- Most notabl y' in Ref - a 

hydrodynamic approach is used to attribute the double 
hump structure in the away side jet to the propagation of 
pressure gradients from dense zones in the plasma which 
can only be observed by averaging over a large number of 
random initial conditions and after subtracting the flow 
component of the azimuthal correlation function. 

However, one can ask whether these features are al- 
ready present at the level of p + p collisions, albeit ob- 
scured by the smallncss of their intensity. For instance, 
two-to-three parton processes would lead to such a dou- 
ble hump structure even in p + p collisions, although 
these would be suppressed with respect to two-to-two 
parton processes. The question is whether in A + A col- 
lisions there exists a mechanism that amplifies hadron 
production from two-to-three processes with respect to 
that from two-to-two parton processes. This question 
has been recently partially investigated in Ref. [14j using 
the event shape analysis in p + p collisions to distinguish 
events containing three jets from those containing two. 

In this Letter, we study angular correlations between 



2 



three hadrons produced both in p + p and A + A colli- 
sions at RHIC using the full Leading Order (LO) matrix 
elements [HI for two-to-three parton processes and com- 
pare these to angular correlations between two hadrons 
using also LO matrix elements for two-to-two parton pro- 
cesses. To the best of our knowledge, this is the first 
study of its kind using the complete LO matrix elements 
to calculate three hadron production in A + A collisions. 
We should mention that the implementation of Next-to- 
Leading Order (NLO) matrix elements for two-to-three 
processes requires the corresponding nuclear parton dis- 
tributions and medium modified fragmentation functions 
to matching accuracy. The former are available in the 
literature [la ], the latter are not fully known. Therefore 



here we use the LO matrix elements, given that we focus 
on the nuclear effects. We use a medium modified frag- 
mentation function in order to treat the medium-induced 
parton energy loss, following the approach of Ref. [HI 
where two-hadron production cross section in A + A col- 
lisions is studied. We argue that since in A + A collisions, 
the final state partons with the short (long) trajectory in 
the away side has a large probability to punch through 
(get absorbed), the effect produces, on the average, a 
double hump in the azimuthal correlation. 

Using the LO matrix elements for two-to-three parton 
processes, including the final state phase space factors 
and enforcing momentum conservation, the three hadron 
production cross section can be written as 
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where yi, ha are the rapidity and magnitude of trans- 
verse momentum of the z-th produced hadron, fa (fa) is 
the azimuthal angle between the second (third) hadron 
and the first hadron, assumed to be the leading hadron 
and y/S is the total center of mass energy. The parton 
distribution functions are denoted by fj/ v (x\, jj?) and are 
given by the CTEQ6 parameterization pjj. u(zi, fj, 2 ) 
are the unmodified fragmentation functions given by the 
KKP parameterization . All along this work, the scale 
p 2 is taken to be the same for all the distribution and 
fragmentation functions and given as the invariant mass 
squared of the final state hadrons, which we take as pi- 
ons. We will focus on the mid-rapidity region and thus, 
hereafter, set all rapidities equal to zero. The momentum 
fractions xi,x%,zi, z 2 are all given in terms of z 3 and the 
transverse momenta of the produced hadrons. Further- 
more, requiring that all momentum fractions arc between 
and 1 introduces the lower and upper limits on the z 3 
integration. Further details will be provided in Ref. [l9l ]. 

We now use Eq. |T|) to calculate the three hadron pro- 
duction cross section in p + p for mid-rapidity at y/S = 
200 GcV. The open circles on the left panel of Fig. [1] 
show the differential cross section as a function of the 
azimuthal angle <fi for the particular case where each of 
the three hadrons carry transverse momentum of h t = 
10 GeV/c and therefore, due to momentum conserva- 
tion on partonic level, they are separated by an angle of 
A(f> = 2n/3 radians. Only the angular distribution of the 
away side hadrons (with respect to the direction of one 
of the hadrons -the would be leading hadron- which in 
the symmetric case chosen here, can be any one of them) 
is shown. 

In order to investigate the effects of a medium on three 



hadron production in A + A collisions, we use the mod- 
ified fragmentation function proposed in [l6| given by 



Dh/i(zi,H 2 



where z[ = 



(1 



x D° h/g (z' g ,n 2 )] +e-^Dl fl ( Zl ^ 2 ), (2) 



is the rcscaled momentum fraction 



-AEi 



of the leading parton with flavor i, z' g = (h)-£jf is the 
rescaled momentum fraction of the radiated gluon, AEi 
is the average radiative parton energy loss and ( 4) is the 
average number of scatterings. The energy loss AEi is 
related to the gluon density of the produced medium via 



,dE, 



dr- 
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p g (r,b_L,r t + nT), (3) 



where b± (note that bt is used to denote the transverse 
momentum of a parton) is the impact parameter of the 
collision. Since here we consider the most central colli- 
sions, we set 6j_ = in the rest of the paper. r t is the 
transverse plane location of the hard scattering where 
the partons are produced and n is the direction in which 
the produced hard parton travels in the medium. The 
average number of scatterings (j), the one dimensional 
energy loss (4§)i<j and the gluon density of the medium 
p g -which is related to the nuclear geometry of the pro- 
duced medium- are taken from Ref. (l6l | where we refer 
the reader also for details on the chosen parameters. 

The three produced partons will in general travel 
through the medium at different angles with respect to 
the direction of r t until they leave the medium after which 
they hadronizc. Therefore, the fragmentation function 
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FIG. 1: Comparison of the differential cross sections for two- 
to-two (full circles) and two-to-three (open circles) processes 
in p + p (left panel) and A + A (right panel) collisions. 



for each produced parton will have a different dependence 
on the path length traveled through the medium. Based 
on simple geometry one can show that out of the three an- 
gles involved, only one is independent which we take to be 
the angle that the direction of motion of particle 1 makes 
with the position vector of the scattering center. The 
normalized cross section is obtained from Eq. (TT]) using 
the modified fragmentation functions given by Eq. 
integrating over the location of the hard scatterings and 
dividing by the nuclear overlap area. Due to the kine- 
matics involved in mid-rapidity collisions, we expect the 
nuclear modification of parton distribution functions to 
be small [2(| HH and therefore ignore them. 

The open circles on the right panel of Fig. [1] show 
the results of the above analysis for the three hadron 
production differential cross section as a function of az- 
imuthal angle </> for mid-rapidity Au + Au collisions at 
\/Snn = 200 GeV. As in the p + p case, we also consider 
that each of the three hadrons carry momentum with 
magnitude h t = 10 GeV/c. Since in our picture frag- 
mentation occurs collinear to the direction of the origi- 
nal parton, the final hadrons are separated by an angle 
Acf) = 27r/3 radians. One can check that the ratio of the 
distributions in the Au + Au to the p + p case is as a 
function of <f> approximately constant ~ 0.7. 

In order to compare the three-hadron to the two- 
hadron production differential cross section in the p + 
p and the A + A cases, we now compute the latter as a 
function of the azimuthal angle (f>. To this end we use the 
LO DIPHOX algorithm 01 . We use CTEQ6 distribution 
functions and in the p + p case, unmodified fragmenta- 
tion functions given by the KPP parametrization. For 
the A + A case, we use the modified fragmentation func- 
tions in Eq. (J5]) integrating over the location of the hard 
scatterings and dividing by the nuclear overlap area. To 
have a fair comparison to the already calculated cross 
section for two-to-three processes, and similar to the ex- 
perimental situation, we take the leading and away side 
particles also having h t = 10 GeV/c. The full circles 
in Fig. [1] show this cross section. On the left panel we 




FIG. 2: (Color online) Distribution of scattering centers (up- 
per panel) and path lengths for two-to-two and two-to-three 
processes (lower panel) for the away side partons, obtained 
by discarding the shortest of the path lengths in each case. 
When three partons are present in the final state, the away 
side parton with the average large (short) path length has a 
greater (smaller) path length than the away side parton in 
the case with two partons in the final state. 



show the results for p + p collisions and on the right 
panel, the results of Au + Au, both for ^/Snn = 200 
GeV collisions at mid-rapidity. Notice that in both cases 
the two-to-three hadron production cross section is sup- 
pressed with respect to the two-to-two result. However, 
and this is the main result of our work, this suppression 
is smaller in A + A collisions than in p + p collisions. 
Dividing this ratio in A + A to that in p + p, we get as 
a function of 4> approximately a constant ~ 2.26. 

Since in the present analysis, the sole ingredient that 
distinguishes between the p + p and the A + A cases 
is the energy loss of partons that hadronize collinearly, 
then the only reason for the cross sections in two-to-three 
processes to be less suppressed in A + A than in p + 
p, when compared to the two-to-two processes, is the 
different geometry for the trajectories of three as opposed 
to two particles in the final state. 

To test this idea, we compute the distribution of path 
lengths with two and three hadrons in the final state. 
We take a nuclear overlap area with a distribution of 
scattering centers denser in the middle and decreasing 
toward the edge. In each case we disregard the shortest 
path length (the one that would correspond to the trig- 
ger particle) and compute the distribution of the other 
path lengths (the ones corresponding to the away side 
particles). This is shown in Fig. O The upper panel 
shows the distribution of scattering centers and the lower 
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one the distribution of path lengths. Notice that indeed, 
when there are three hadrons in the final state, the large 
(short) path length in the away side is greater (smaller) 
than the case of the away side particle when there are 
two hadrons in the final state. This means that for the 
case of three particles in the final state, even if one of the 
non-leading particles gets absorbed by the medium (the 
one with the largest path length) , the remaining particle 
has a large probability of punching through than in the 
case when one has two particles in the final state. We 
interpret this as signaling that in two-to-three processes, 
there is a large probability to have one of the two away 
side particles being absorbed an the other randomly get- 
ting out, producing on the average, a double hump. 

In conclusion, we have presented a LO calculation in 
p + p collisions for the cross section of two-to-three pro- 
cesses and compared to the corresponding cross sections 
for two-to-two processes, using the same approximations. 
We have also computed the same processes in A + A col- 
lisions taking into account a detailed Glauber approach, 
similar to the one employed in Ref. [23[ . The results show 
that the cross section in two-to-three processes in A + A 
is less suppressed -with respect to two-to-two processes- 
than in the case of p + p. We attribute this result to a 
purely geometrical effect associated with the differences 
in the path lengths of three as opposed to two particles in 
the final state in the case of A + A collisions. These re- 
sults raise several interesting possibilities that need to be 



addressed, among them we can point out the following: 
1) Two-to-three processes are naturally expected and the 
strenght of their signal is not beyond the capabilities of 
current experiments. 2) Considering that two-to-three 
processes exist in p + p collisions, as suggested also in 
Ref. [Til , our work shows that their observation in A + A 
should be enhanced with respect to the strongly energy- 
loss suppressed two-to-two processes and thus, that this 
effect may have a bearing on the shape of the away side 
for different kinematical cuts in azimuthal correlations in 
A + A. A detailed study of this effect is currently under 
way and will be reported elsewhere. 
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